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ABSTRACT 


A numerical method for two-dimensional cavitating flow 
is developed for the flat plate. The linearized boundary 
value problem is restated as a set of coupled integral equations. 
The integral equations are approximated numerically. The 
numerical approximation is executed by a Fortran IV computer 
programe The computed results are compared to the analytic 
solution. This method should provide insight into developing a 
method for three-dimensional cavitating flows and is readily 
extendable to cambered profiles. 
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NOMENCLATURE 


Cavity Area 

chord length 

coefficient matrix 

element of coefficient matrix 
pressure coefficient 

lift coefficient 

location of vortex control points 
location of source control points 
locus of points describing mean camber line 
and cavity surface 

cavity length 

number of vortex elements 

number of source elements 
pressure 

pressure at infinity 

velocity 


source distribution 


= perturbation velocity components 


inflow velocity 

complex velocity = u-iv 

coordinate system 

lower boundary of vortex or source element 


complex plane = x + iy 





a 
( ) 
io) 
o)- 


Dy 


XG) 
§ 7 


s+ 8 COA Dd 


quantity with a bar is a dimensional quantity 
no bar indicates a nondimensional variable 
variable to be evaluated at x = 0* 
variable to be evaluated at x = 0” 
variable to be differentiated wrt x 

angle of attack 

vortex distribution 

mapped plane = a +0 7 

dummy coordinate system for integration 
and see C above 

fluid density 

cavitation number 

velocity potential 


perturbation velocity potential 


= degree of singularity 





CHAPTER I 
INTRODUCTION 

The objective of this investigation is to develop a 
numerical solution to cavitating flows which can be compared to 
a known analytic solution. Therefore, the cavitating two-dimen- 
sional flat plate in steady flow has been chosen for this 
investigation. Geurst [1,2,3] has given solutions for the 
linearized problem for both partial and super cavitation. 

The purpose in considering this problem is to determine 
and overcome the numerical difficulties associated with this 
problem. Adequately solving the two-dimensional problem is a 
preliminary step in developing a numerical solution for the far 
more difficult three-dimensional and unsteady flow problems. 

Further the numerical method can be applied to foils 
of different camber lines. Geurst [2] has obtained a camber 
line solution for the case of parabolic camber in partial 
cavitation only. 

Widnall [4] developed a numerical method for three- 
dimensional, unsteady, supercavitating flows. However, Widnall 
used assumed cavity lengths and demonstrated that, for long 

cavities, lift coefficients were relatively insensitive to 
| cavity length. Clearly for short super cavities, on the order 
of chord length, and partial cavities, a proper statement of 
cavity termination is required for the proper solution of 


cavitating flows. 
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In this respect this report is a preliminary investiga- 
tion into the inclusion of cavity termination in cavitating 
flows. The results of this report should be applicable to 
three-dimensional and unsteady cavitating flows. 

The solution of two-dimensional, cavitating flows pre- 
sented herein envolves five distinct steps. These are: 

1) statements of the linearized boundary 


value problem, 
2) solution of the boundary value problem in 
the form of coupled integral equations in 


source and vortex distributions, 


3) a numerical approximation to the integral 


equations, 


4) execution of the numerical approximation via 


FORTRAN computer progran, 


and 


5) validation of the numerical results. 


11 

Linearizing the boundary value problem and rewriting 
the problem in the form of integral equations is not the only 
method of solution. Another possible approach is the finite 
difference method. Jeppson [546 and Mogel and Street [7] have 
taken the finite difference method as far as a circular disk 
perpendicular to the main stream with an infinite cavity. The 
major disadvantage of the finite difference method is the large 
programing and computational effort required. 

The integral equation method for the linearized problem 
represents a far smaller computational effort. It is for this 
reason, plus the existence of the analytic solution, that the 
integral equation method is used. For an excellent discussion 
of the various methods of solution see Birkhoff (3) e 

The following chapters of this report will discuss the 


stated steps of the solution. 
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CHAPTER II 


THE LINEARIZED BOUNDARY VALUE PROBLEM 


Geurst {1,2,3} has presented a development of the 
linearized boundary value problem for both partial and super 
cavitation. A similar development is given here, except that 
the coordinate system and the nondimensionalization of variables 
differ. These differences are strictly a matter of convenience. 
Geurst*’s method of solution relies on a series expansion of the 
complex velocity in the neighborhood of infinity. This in 
turn requires that the boundary value relations to local slope 
on the foil and cavitation number be stated at infinity. 

The neighborhood of infinity is not convenient to 
computer solutions. Therefore the boundary value relations to 
local slope on the foil and cavitation number need to be stated 
on the foil and cavity. 

The coordinate system and the foil-cavity relation to 
the coordinate system are shown in figure 1(a). The inflow 
velocity, U, , is assumed to be parallel to the x axis and 
uniform in the y and z directions. The angle of attack, oh , 
is taken to be the nose-tail line for cambered profiles. 

First the steady flow form of Bernoulli’s equation is 


written between infinity and another point in the flows 


2 Piist g=elvd 
(1) hobaeeat’, 371" 9/ 


Dil gd! 
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This equation is now rewritten in the form of a nondimensional 


pressure coefficient. 





B-P gs 
72 = ( os gal) 


(2) = 
Cp + 


“vO! 


On the cavity the pressure coefficient is the cavitation 
number, (GJ »« The velocity potential can be written as the sum of 


the potential due to the inflow velocity plus a perturbation 


potentials 


(3) d=U,X +P 


Where 


——_—o 


D = perturbation velocity potential. 


From equation (3) the x and ¥ components of velocity are: 


0M ~,WPi_oea 
(4) ab 60,4 ae 


ag 


C5) 


ww) QV 
RA 
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Where 


ol 
ti 
4 


x component of the perturbation velocity 


<\ 
fl 
< 


y component of the perturbation velocity. 


It is assumed that the perturbation velocity is very much less 


than the inflow velocity, that is: 
UU << Uso 


On the cavity we now have: 
ek — = = —— 
Q0=[02 +240,+42 + G7 
using(6) 
me 
Be _ ee 
Thus equation (2) for the cavity becomes: 


(8) ud. = T/ 2 


Where 


and (4, is now used as a nondimensionalization constant. 
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Equation (8) is the statement of the linearized boundary con- 


dition on the cavity surface. The next step is to obtain the 


boundary condition on the wetted surface of the foil. 


On the foil wetted surface the boundary condition is 


flow tangency. 


vr e.da 
(9) Unjt4 gx 


u aT | _ ave 
Bom ofa o. fe higher order terms = T¥ 
or 
ao 3 dh 
(/0) Se a 


Also, it is assumed that the foil is very thin compared 
to the chord length such that the boundary condition in equa- 
tion (10) can be applied to the mean camber line. 

Further the entire foil cavity system is collapsed to 
the X-Z plane and the boundary conditions applied on the x-2 
plane with all length dimensions nondimensionalized by the chord 
length, figure 1(b). Another condition to be satisfied is cavity 
termination. 

Geurst [2] proved that for the linearized problem, the 


re-entrant jet and Riabouchinsky models for cavity termination 


reduce to a statement that the cavity closes at its end. This 


condition can be written as 


2 16 
» fle 
M) I Cavity -( iz a | =0 an 


or 9 
(12) [ LY — h, | = 


The next condition to be satisfied is the condition at 
the trailing edge of the foil (x = 1). This condition requires 


that there be no pressure jump accross the foil at the trailing 


edgee 
13 Ce) —C500) =O 


In summary the conditions that the solution must 


satisfy are: 


C/ r- 
se TL 2 OS 2 


» on wetted surface 
(C2) ye ns of the foil 


(C3) [Tx - | =0 


(C4) Cp (1) GG) =O 
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The last condition is that the perturbation potential 


satisfies Laplace's equation: 


(C5) RP ease Bn 


The next chapter will give a solution to Laplace's 


equation which satisfies the above boundary conditions. 
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CHAPTER III 


THE INTEGRAL EQUATION FORMULATION 


This chapter presents the development of the integral 
equation formulation of the problem. This development proceeds 
from the solution for the velocity field induced by a two- 
dimensional distribution of vorticies and sources in space. 
Then it can be shown that this velocity field will satisfy the 
boundary conditions. This results in a set of coupled integral 
equations. 


First consider the flow to be in the complex (2) plane, 
Where 3 


Z2xX+ ly 


and the complex velocity is given by: 
ui(Z)2 U—-liv 


Then the appropriate general form of the solution is [9] 


fae «ur(Zz) = eas —— de — ae 2 
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or 
| a | 
(15) weay=an | “EEGD de 


where, 


| G7 O *) ie © 


9 (E)= source distribution over the cavity 


YCE)* vortex distribution over the foil 


1 for a partial cavity 
+ for a super cavity 
g¢§)= 0 for £<§ for a partial cavity 
and 


a(§) = 0 for 1<§ for a super cavity. 


Plemelj*s formulas can be applied to equation (15) 
as follows [10]. 


(16) WTC) — W(x) = - YX) — 63902) 


a 
— ¥(0€)-¢3C§) 
17) Wtix) + were spp GOS dg 
oO 
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When equations (16) and (17) are added together and 


the real part is taken, the result is: 


XB 
BOF) 


oO 


or, with the cavity boundary condition, 


(/8) Ga = f(x) t +f $d. wa oS 


with 
W*AV=z=O foy x2] 


The integral in equation (18) is a Cauchy principal value in- 
tegral over the length of the cavity. 
When equations (16) and (17) are subtracted from each 


other and the imaginary part is taken the result is: 


/ 
ca (x) = > Be dat a d¢ 
O 


or, with the wetted surface boundary condition, 


/ 
= Cz) / X¥CE) 
(19) h(a - BF 4 anf XD as 
with ° 


at = OC for X> 
Again the integral in equation (19) is a Cauchy principal 
value integral. For a flat plate at an angle of attack 


(20) h( X)=x—-aX 
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and this is the case taken throughout the rest of this report. 


The source distribution q(x) represents the slope difference 
between the cavity surface and the mean camber line of the 
foil at x. While the vortex distribution Y (x) represents the 
difference in the x component of perturbation velocity between 


the upper and lower surface at xX. 
ee 20%) = V(X) — V (2%) 


eee (4%) = — (K+ U(X) 


These statements follow directly from equations (16) and (17). 
When equations (2), (7) and (22) are combined the jump in 


pressure on the foil iss 


G3) Cp (®) - Cy (x) = -a YX) 


and the coefficient is: 
/ 
aoe Cc, - -a | Yooda 
O 


The combination of condition (C 3) and equation (27) gives, for 


the closure condition, 


2 
(25) [aca 26 


O 
It is clear from equation (23) that condition (C 4) is: 


@6) Y¥(1)=0 
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Equations (18), (19), (20), (25) and (26) constitute a complete 


problem statement for the steady flow cavitation of a two- 
dimensional flat plate, only, it is convenient to rearrange 
equations (18) and (19) into a form that contains the ratio 


of angle of attack to cavitation number, 


Q7 4=- *@ 4 oft acevo de 


LB = Opole BO) YE 
@8) Oo anf Le oY 


Equations (27), (28) and (29), subject to Y(=oO » are 
sufficient to obtain cavity length, the source distribution 


and the vortex distribution as functions of the ratiox/c only. 
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The next section will develop the numerical approxima- 


tion to equations (27), (28) and (29) and the method of solution, 
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CHAPTER IV 
NUMERICAL APPROXIMATION TO THE INTEGRAL EQUATIONS 


The method used in the numerical approximation is to 
approximate the integral equations (27), (28) and (29) by a 
finite set of linear algebraic equations with N unknown source 
densities, M unknown vortex densities and the ratio «</o ° 
The numerical integration procedure requires that cavity length 
be the independent variable and = be computed. 

This method requires an assumption about the functional 
form of the source and vortex distributions. It is assumed that 
the vortex and source distributions are constant over a small 
element of the foil or cavity, figure 2. The last vortex 
element is assumed to be linear with a value of zero at the 
trailing edge to satisfy equation (26). The amplitude of each 

% » G is unknown. 

In each vortex element a control point is placed at 
which the flow tangency boundary condition is satisfied. In 
each source element a control point is placed at which the 
cavity boundary condition is satisfied. The vortex control points 
are at X; and the source control points are at 2, . The 
relations for the control point locations are: 


1) for the vortex control 


(30) Cid 72; r ( Xeje, i Aaj) F 5 DET 
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and 
2) for the sources 
(BI ys = Xyet (Xin — XIE o<f <I 


The fractions f and ff are yet to be determined on the 
basis of obtaining stable results. 

There will be one equation for each vortex control 
point, giving M equations. There will be one equation for 
each source control point, giving N equations. Plus there will 
be one equation for the cavity closure condition. The vortex 


equations are: 


Xe jt 
Ga) ofan —vtiyf “te ye 4 


Xie) / 
— & S/T ot > 
) Lae Xj-§& dg ela OG ~ ELI — Xpy) d§ 
Xen 


and for the a control point 


as 
rae “tt (1- ¥) 
(3) O=an=S mea Cy oa (Xp EC - Xey) dé fe 
Kip4 


[Se y. 
7 Ou d 
5 
Where, 


ee for J7N and Bf, =0 for MIN 


| 26 
The source equations are: 


N Abie 
os i 
G4) W=-E& +) [ee de + 
Tia 4 





x, “I~ 9 G=O for gom 


and the closure condition 
AQi+i 
Be 
ie) ) L ad O 


Equations (32) through (35) are the numerical approximation to 
the integral equations. When the indicated integration is per- 
formed equations (32), (33), (34) and (35) become 
M-l 
$5. YX ita Yi; ee 
— a7 — — + —L ——— a xj Le 
iz} 


+ St | po Ln ( Bae aa +/J 3 $= 0 for j>v 


37) o=an& ~p B44 Ym [y+ Ci-F) Lo (Fe)] + 


M-l 


Xm ~ Xe 
Ln (F Xgivs) 5 $4 =0 DT A ee hd 


/ 


a 


(38) 


and 


(39) 
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ee = 
. 7 + T Be bn (ee x) yg = bn (5H) 
3} 


Xi O ee J>N 


A/ 
» aoe (Kein = X92) = O 


L=/ 


Equations (36) through (39) are linear algebraic equations in 


ee zs It is convenient to rewrite 


OX 
the unknowns Go an e 


these equations in matrix form. 


(40) 


oi O 
or 
ne 
oT O 
4. 1) = 

|=)" 
9 
cf Tt 
& O 
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The coefficient matrix,(C], is formed from the coeffi- 


cients of = , ee 9 as . a and ui in equations (36) 
through (39). 
The first submatrix, (I), is formed as follows by the 


coefficients of % in equations (36) and (37). 


_ Zz 
(4) Cj, = Ln (TF , JtM 


(42) Cum= 1+ U-#) Ln (4, 


Xe¢ : 
(43) Ci 7 = bn ( ae ae 5) CHEM and j#mM 


mee xj — Xe 
(44) Cim= Foy Ln( ZS) + , 





J#I™ 


The second submatrix, (II), is given by the coefficient 


of 9j/¢ in equations (36) and (37): 


(45) Ciyjem=-W, feign 


(46) Ci - 0, aA Mtlétl ¢Mrn 


and L£eytm., 
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Submatrix (III) is the coefficient of «~/7 in equations 
(36) and (37): 


(+7) Cj mrt eds for /<j<M™ 


(48) Cj) ment =O for Mey 4 MtNt/ , 


The source equations are contained submatricies (IV), 
(V), and (VI). The submatrix (IV) contains the coefficients 


of %% in equation (38) : 
and 


(50) € = 0 otherwise. 


Submatrix (V) contains the coefficients of & in 


equation (38): 


‘ Xj — Kec : : 
 Cojmec = Pomme) 9 Mtl ete 


LSJ teNn 


and 


(52) Cmtjmt = = Ln (te —f{; ) 


The last 
equation (39): 


(53) Cueneyi aia XOcrl ae Ker 5 Mtl = l = MtN 


Lse 
Cpanel é O otherw 
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submatrix (VI), contains the closure condition, 


To solve equation (40) for the unknown vortex and source 


distributions and the ratio ~¢@ the coefficient matrix is ine 


verted. 
the boundary condition matrix. 


Thus the solution can be written ass 


¥, oO 
¢ 
al _ ZI es 

(64) fal — C rr 
By 


0 2: 


The inverted coefficient matrix is then multiplied by 


31 
On the surface this procedure appears to be straight- 


forward and should produce a solution without difficulty. 
However, this is not the case. 

Difficulties in obtaining a meaningful numerical solution 
lie in the selection of control point locations within each 
source or vortex element. Control point locations are critical 
to approximating the singular behavior of the analytic solution. 
It was found that vortex control points should be at 90 percent 
of the element length and source control points at the element 
midpoint. However, special care must be taken with a few of the 
source control points. 

The critical control points are the first and last 
source control points. A systematic variation of the control 
point position determined that the last source control point 
should be placed at 10 percent of the element length. As an 
example of the importance of this control point, the variation 
of the computed ratio of = with control point location is 
shown in figure 3 for a cavity 0.5 times chord length. 

The selected control point position is at 10 percent 
of the element length. Variation of the first source control 
point indicated the placement should be at 90 percent of the 
element length. Since the source distribution is singular at 
the cavity leading edge and termination, these control points 


should be placed away from the singularities. 
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Variation of other control points (second and next to 


last source, last vortex, first source past trailing edge for 
super cavity) showed the best placement to be the generalized 


locations. Thus the fractions f and f, ares 


f = 0.90 for all vortex elements, 

f, = 0.90 for first source control point, 

fg = 0.10 for last source control point, 
and 

fg = 0.50 for all other source control points. 

The computer program that executes the solution to the 
described numerical procedure is listed in Appendix A. This 
program is composed of a main program, which performs Input/ 
Output operations and logic control, and three subroutines. 

The first subroutine, CPGEN, determines the vortex and 
source element sizes and the locations of the control points. 
The required input information for this subroutine isthe cavity 
length, the number of vortex elements and the number of source 
elements. 

The second subroutine, MATRIX, computes the coefficient 
matrix, LC], and the boundary condition vector. The required 
input information for MATRIX is the output of CPGEN and the 
number of vortex and source elements. 

The last subroutine, RMINV, inverts the coefficient 
matrix. This subroutine is a standard MIT matrix inversion 


routine({11]. The inverted matrix is then multiplied by the 


3 
boundary condition vector, in the main program, to obtain P 


od VE Be 
the solution vector for = ° = and = . 
A special note on the program notation is made here, 
The vector XU(I) is the upper boundary of a vortex or source 


element. This vector is given by: 
XU(I) = XL(I + 1) 


through an assignment statement in CPGEN. 


CHAPTER V 3 


ANALYTIC SOLUTION 


The method of gauging the numerical results is the 
analytic solution for the cavitating flow of two-dimensional 
flat plate. The first half of this chapter will discuss Geurst's 
Solution for partial cavitation [1,2]. The second half will 
discuss Geurst’s solution for super cavitation [3] . 

Geurst solves the partial cavitation problem by a 
conformal mapping of the 7 plane to the C plane as follows. 
First the 2% plane is given bys 


+ corresponds to leading edge 


See eo = 
Gs) Zz ot > Ze +! correspends to traclingedge 


The mapping function iss 


eects 5 GHS HR 


This mapping is such that the point at infinity in the physical 
plane is mapped to (=i » The flow is then in the upper half of 
the ¢ plane. Figure 4(a) shows the resulting plane. 

The solution for the complex velocity is: 


(57) ws) _ 2¢+é 


oT 
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where, = 
my lee 
i 7 
= 4 eog— 
, oo 
Bz [tang Hee sinc F-£)- cos(% - $)]l aime] 
and 
= [ton 5 toting Cost —£) + sin FS aver] 


With equation (57) is substituted into equations (21) 
and (22) over the cavity and wetted surface of the foil, the 
vortex and source distributions can be obtained. Figures 5 
and 6 show the vortex and source distributions for cavity equal 
to one-half chord length. These figures show the leading edge 
and cavity termination singularities. One of the difficulties 
for the numerical procedure is to reasonably approximate these 
four singularities. 

Other quantities of interest are the ratio = » the 
lift coefficient and the cavity area. These are given as: 


J-stin§ 


68) La btans (Ms 


aoe 
= 
C 





Pe ona lit nF | 


ahd 
(60) Ares — Beg Fier s2e08 cosd sind + 





271K 
f + cotand(l+sin$§)(it 35th $§—-Aasin2 §— é sing) 


36 
The equation for cavity area given in/{1] is in error 


and was corrected to equation (60) in [2]. Plots of equations 
(58) through (60) are shown in figures 7, 8 and 9. 
To obtain a solution to the super cavitating problem, 


Geurst performs a similar mapping (figure 4(b) )defined bys 





(6/) ee q-% 


The solution for the complex velocity is then. 


iS 
(62) wet) girs \/ * 








where 
b= Gop = cotan’d 
+= COS 5 
at. v5ind [ Zang sinC#-£) — cos(%-$) | 
o aX 
and 


Beeisins [ tan § cos(#-#) + Sin (F -3)/ 
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The vortex and source distributions are obtained using 
the same procedure as with the partial cavitation case. Figures 
10 and 11 depict these distributions. In the super cavitation 
case only a leading edge singularity appears in the vortex 
distribution. Again the source distribution shows the leading 
edge and cavity termination singularities. 


The quantities -“ , lift coefficient and cavity area 





are: 

ee pag 

C, / 

SS OEE 

(64) IM a asin&(lt sinS) 
and 

Area _ Mees. 
(65) ATTY X 3a sin §(i—sin8)* 
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Figures 12, 13 and 14 depict equations (63) through (65). 


As with partial cavitation when the cavity approaches chord 
length the lift coefficient and cavity area become singular. 
In chapter VI the computed results are compared to the 


analytic solution. 
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CHAPTER VI 


RESULTS, CONCLUSIONS AND RECOMMENDATIONS 


The primary result of this investigation is a numerical 
method for two-dimensional cavitating flows which gives computed 
values close to the analytic solution. The computed results 
have been plotted in figures 3 and 5 through 16. In figures 5 
through 14, the analytic results of Geurst are also shown. 

Figures 5 and 6 show the computed vortex and source 
distributions for a cavity of one-half chord length. Each com- 
puted point is plotted at the control point location, except 
for the last vortex element which is plotted at x»,,. It is to 
be remembered that each computed point represents a constant 
value from Xg, to X,,,,. These two figures show that the com- 
puted distributions fit the analytic solution reasonably well. 
However, the largest magnitude elements are not shown on the 
figures (see the sample computer output in Appendix B for their 
values). These elements do not appear to approximate the 
proper singular behavior. 

Near the singularities the distributions of sources and 
vorticies can be approximated by 1/x", For the leading edge 
n has the value one-quarter and for cavity termination n 
has the value one-half{1]. The computed values of n are 
given in Table 4 for three values of M and Ne From Table 4 it 
iS apparent that the behavior of the two closest elements does 


not match closely the analytic solution's singular behavior. 





4oO 
Figures 7, 8 and 9 compare the computed values of «“/7 , 


Cjaypn and A/ypy, » for ten cavity lengths to the analytic 
solution. The plotted results show the points closest to the 
analytic solution from a number of computations for various 
values of M and N. In general the computed results compare 
quite favorably with the analytic solution. The exception to 
this is for cavity lengths approaching chord length. For these 
cavities the method appears to be very slowly convergent. For 
these cavities a large number of elements are required for the 
cavity-foil region leading to a large number of equations. 

Table 3 contains all the computed results to date. This 
table gives results for variations of M and N for specific 
cavity lengths. For example, taking £= 0.5, the results show 
that when the vortex and source element size is smaller in the 
combined cavity-foil region than in the fully wetted region, 
the computed values show an improved level of convergence to 
the analytic solution. This result is, also shown in figures 
5, 6, and 15. . 

Figure 15 Shows the results of a convergence test for a 
one-half chord length cavity . this figure also shows the 
tendancy for better convergence for a smaller interior element 
size. It is noted that for figure 15 the vortex element size in 
the region of fully wetted flow are the same for both ratios 
of M and N. This effect is also shown in figure 5 for the vortex 
distribution with the same ratios of M to N and the source 


distribution in figure 6. 





This conclusion points out one of the problems with 
the numerical method used herein. The method reported on requires 
that the vortex and source elements in the combined cavity- 
foil region be of equal size. An independent variation of the 
source and vortex element sizes cannot be performed. Thus it is 
not possible to determine whether decreasing the size of the 
vortex elements, the source elements,or both, results in the 
improved convergence. It is recommended that an investigation 
of effects of independent variations of the vortex and source’ 
elements be done. This can be achieved through a moderate re- 
vision of the computer program. 

Also, since the singularities do not appear to be well 
matched, a further investigation of the location of control 
points should be considered. It seems desirable to have an 
analytic basis for the control point locations. 

From the general form of the computed vortex and source 
distributions, this investigator is convinced that assuming 
piecewise linear distributions will produce an immediate improve<- 
ment in the computational results. It is, therefore, recommended 
that this numerical method be modified to incorporate piece- 
wise linear vortex and source distributions. 

The results for the super cavitation case are very 
Similar to the partial cavitation case. The convergence to the 
analytic solution becomes worse as the cavity length approach 
chord length. The same conclusions and recommendations made 


for the partial cavity apply to the super cavity. 
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Another purpose of this investigation has been to 
develop a method of allowing arbitrary values of angle of attack 
and camber as inputs and then determine the cavity length. An’ 
itterative procedure to accomplish this has been developed. 

The method and the computer program are in Appendix C with the 
computed results. The method and program are a preliminary effort 
provided to prove only that such a method is possible. 

In summary these are five recommendations for future 
investigation. 

1) <A program be written that allows the vortex 

and source element sizes to be determined 


independently. 


2) An analytic and further numerical investigation 


of control point locations. 


3) Piecewise linear distributions of sources 


and vorticies should be investigated. 
4) That the method herein not be applied to the 


three-dimensional case without further 


investigation. 
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5) 


6) 
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The program should be modified to include 


cambered profiles. In this case the ratio 
of “/7g cannot be computed as an unknown. 
Either the angle of attack or cavitation 
number must be known along with the camber 
line. It is not clear that the control 
point locations will be the same as with 


the flat plate. 


Further development of the method of 
computing cavity length from arbitrary 


values of ro and om e 


The conclusions drawn from this investigation ares 


1) 


2) 


3) 


the method gives computed values close to 


the analytic solution, 


the distribution of sources and vorticies 


is probably the best test of computed results, 


the method shows poor convergence for 


cavity lengths near chord length, 





4) an iterative procedure to determine cavity 
length from an input of angle of attack 


and cavitation number is possible. 
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Table 1 
Values of M and N for figures 5 through 9 


No. of Vortex No, of Source 
Cavity length elements 
Q M 


0.1 
0.2 
0.3 
04 
0.5 
0.6 
0.7 
0.7 
0.8 
0.9 
0.9 
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Table 2 
Values of M and N for figures 10 through 14 


No. of Vortex 
Cavity length elements 


60 








61 
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number of equations 
figure 15 
Numerical Convergence, £=0.5 
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cpu time 
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50 40 60 80 
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figure /6 
Computation time 





Table 3 


Computed Results 





Points plotted in the figures. 
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elapsed 
time 


(pp $8¢) 


4. 66x1072 
3 »45x1072 
.0129 
60119 
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Table 3 continued 
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Table 4 


Singular Behavior 






n* 


termination | 
source 








termination 
Vortex 





* Based on only the two elements closest to 
the singularity (e.g. for the leading edge 
vortex distribution the elements are the 
first to vortex densities) and @«= 0.5. 
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Appendix A includes the flow charts for the computer 
program and the computer program listing. The flow charts 
are provided to give a broad overview of the program flow. 
The numbers between blocks in the flow charts indicate 
approximate statement numbers in the program listing. 
A flow chart for RMINV is not provided since it is a standard 


matrix inversion routine, 
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APPENDIX B 
SAMPLE COMPUTER OUTPUT 
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Appendix B contains the sample computer output. This 


output contains the results of a convergence test for 
cavity lengths equal to # chord length and 1.5 times chord 
length. Note that for a partial cavity source control 
points beyond cavity termination are listed in the output. 
These control points are superfluous and not used in the 
program. The same is true of control points for the 
vortex elements beyond chord length in super cavitation. 
The listing of source and vortex densities has the 


same order as the listing of control points. 





NO. CF VORTICIES= 10 
NO. OF SOURCES= 5 
CAVITY LENGTH= 0.50009 


**¥e DETERMINANT=-.1317091E+09*** 


LIFT COBPICIENT/ (2*PI*ALFA) =0.1092649F+01 
AREA OF CAVITY/(2*PI*ALPA) =C.4117717E-01 
TOTAL SOURCE STRENGTH/SSIGMA =C.2272427E-06 
ANGLE OF ATTACK/SIGMA =0. 1067560E+00 
CAVITY LENGTH =0.5000000E+00 


ELAPSEC TIME= 2 1/100 SEC. 
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XL(T) 
0.0 
Oot 2 505-01 
0.62500E-01 


NO. OF VORTICIES= 20 
NO. CP SOURCES= Te 
CAVITY LENGTH= 0.50000 


X (I) 
0.28125E-01 
0.59375E-01 
0. 11875E+00 


XU (I) 
0.31250R-01 
0.62500E-01 
C.12500E+00 


0.12500£+00 0.18125E£+00 0.18750E+00 
0. 18750E+00 0.24375E+#00 C.25000E+00 
C.~25000E +00 0.30625E+006 0.31250E+00 
0. 31250E+00 0.36875E+00 C.37500E+00 
0.37500E+00 0.43125E+00 0.43750E+00 
0.43750E+00 0.46562 E+00 0.468758+00 
0.46875E+#00 0.49687E+00 C.50000E+00 
0.50000EF+00 0.52812E+00 0.53125E+00 
02531255+00 0.5593 7E+00 C.56250E+00 
0.56250E+00 0.61875F+00 C.62500E+00 
0.6250CE+#00 0.68125E+00 0.68750E+00 
0. 68750F+C0 0.743755F+00 C.75000F+00 
0.75000F+00 0.80625EF +00 0.81250E+00 
0.81250£E+00 0.86875 £+00 0.87500E+00 
0.87500E+00 0.93125E+00 0.93750E+00 
0.93750E+00 0.96562 E+00 0.96875E+00 
0.96875E+00 0. 99687E+00 C.10000E+401 
*** DETERMINANT=0. 35329515+17%** 
VORTEX DENSITIES/SIGMA: 
-0.93763E+00 -O.87444E +00 -0.79342E+00 
-0.6856CE+00 ~0.66212E+00 ~0.63783E+00 
-0.64077E+00 ~0.55656E+00 0.65990E+00 


=U 2997 1E-0'T 
-0.51625E-01 


-0.63426E-01 


-~C.37786E-01 


SOURCE DENSITIES/SIGMA: 


0.80515E+00 0.53662E+00 C.33202 E+00 
-0.5600 7JE+00 -~C.13084E +01 


LIFT COEFICIENT/ (2*PI*ALFA) =0.1137383E+01 


AREA CF CAVITY/s(2*PI*ALFA) =C.4946671E-01 
TOTAL SOURCE STRENGTH/SIGMA =0.1750886E-06 
ANGLE OF ATTACK/SIGMA =C0.9841657E-01 
CAVITY LENGTH =€.5000000E+00 


PesPSsetD TIME= 15 1/100 SEC. 
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XS (T) 
0.28125E-01 
0.46875E-01 
0.93750E-01 
0. 15625E#00 
0.21875E+00 
0. 28125E+00 
0.343 75E#00 
0. 40625E+00 
0.45313E400 
0.47187E+00 
0. 5156 3E#00 
0.54688E+4+00 
0. 59375E+00 
0.65625E+00 
0.71875 E+00 
0. 78125E+00 
0.84375E+00 
0.90625F+00 
0.95313E+4+00 
0.98438 F+00 


-0. 724 74E+00 
~0.65261E+00 


0.11008E+00 


-0.70911E-01 
-0.38749E-01 


0.20736E+00 


-0.30393E+00 





OF VORTICIES= 


OF SOURCES= 


CAVITY LENGTH= 


XL (TI) 
oe) 
0.62500E-01 
Ge 125N0E+00 
0. 25000E+00 
0.37500E+#+00 
0.50000E+00 
0.62500E+00 
C.7500CE+00 
0. 875008+00 
0.937508F+00 
0.10000E+01 
9.10313E+01 
0.10625E+01 
0.11250F+061 
0. 1187556+01 
C.125008+01 
0. 13125E+01 
0.13750E+01 
0.14375E+01 
0. 14688E+01 


10 
29 
1.50000 


X (T) 
0.56250E-01 
0.11875E +00 
0.237508 +00 
9. 36250E+00 
0.48750E+00 
0.61250E+00 
0. 73750£+00 
0.86250E+00 
0.93125E+00 
0.99375E+00 
0.10281E+01 
0. 10594E+01 
0.11187E+01 
0.11812£+01 
0.12437E+01 
0.13062E+01 
0.136872+01 
0.14312E+01 
0.14656E+01 
0. 14969Et+01 


XU (T) 
0.62500E-01 
0. 12500E+00 
0.25000F+00 
C.37500E+00 
0.50000E+4+00 
0.62500E+00 
C.75000E+#+00 
C.87500E+#00 
¢€¢.93750E+00 
0.10000F+01 
0.10313E+01 
C.10625F+01 
0.171250E+01 
0.11875 F+01 
0.12500E+01 
0.13125E+01 
C.13750E+01 
0. 14375E+01 
0.14688E+01 
C.15000E+01 


**¥*k DETERMINANT=0.38512468+13*** 


VORTEX DENSITIES/SSIGSA: 


~0.11499F#01 -0.10381E+01 -C.87268F+00 
~0.62812E+00  -0.54372E+#00 -0.46015E+#00. 
-0.2608&8F+C0 —-C.25544 £400 
SOURCE DENSITIES/SIGMA: 
0.132678+01 0.96 160£+00 0.73057EF+00 
0.50964E+00 C.43215E+00 C.36176E+00 
0. 25084EF+00 0. 23408E+00 0.19550E-04 
-0.16451824+00 -6.37803E+00 -0.49625E+00 
-0.98532E+00 <-0.15239E+01 -6.29826E+01 
LIFT COEFICIENT/ (2*PI*ALFA) =0.5062250E+00 
AREA OF CAVITY/(2*PI¥ALFA) =C.2296074E#00 


TOTAL SCURCE STRENGTH/SIGMA =-.1788139E-06 
ANGLE OF ATTACK/SIG MHA =0.3843980E+00 
CAVITY LENGTH =0. 1503000E+01 


ELAPSED TIME= 15 1/100 SEC. 
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XS (1) 
0.56250E-01 
0.93750E-01 
0.18750E+00 
0.31250EF+00 
0.43750E+00 
0.56250E+00 
0.68750E+00 
0.81250E+00 
0.90625E+00 
0. 96875E+00 
0. 10156E+01 
0.10469E+01 
0.10938E+01 
0.11563E+01 
0.12188E+01 
0. 12813E+01 
0.13438E+01 
0. 14063E+01 
0. 14531E+01 
0. 14719E+01 


-0.72748E+00 
-0. 3701 4E+00 


0.60112E+00 
0.28510E+00 
~0.33281E-01 
-0. 75565E+00 
-0.56086E*01 





NC. OF VORTICIES= 15 
NO. OF SOURCES= 30 
CAVITY LENGTH= 1.50000 


*x* DETERMINANT =0.7520241E+19*** 


LIFT COEFICIENT/ (2*PI¥*¥ALFA) =0.5145274E+00 
AREA OF CAVITY/(2*PI*ALPA) =C.2330343E+00 
TOTAL SOURCE STRENGTH/SIGMA =C0.1430511E-05 
ANGLE OF ATTACK/SIGMA =0. 37883252400 
CAVITY LENGTH =0. 15000008401 


ELAPSEC TIMNE= 39 1/100 SEC. 
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APPENDIX C 
COMPUTER PROGRAM FOR COMPUTATION 
OF 
CAVITY LENGTH FROM ARBITRARY VALUES 
OF 
ANGLE OF ATTACK AND CAVITATION NUMBER 
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The iterative method for computing cavity length from 


arbitrary values of angle of attack and cavitation number is 
based on assuming an initial cavity length ( 2 = 0.5 ) then 


computing the next cavity length as: 


= Pinpot” . 
(IC) ee = BY, (S): 
GT 4¢ 
In equation (1 C) the value of (&), is the value computed 
from {; . This works well for moderate length partial cavities. 
For super cavities the method used is to set upper 
and lower boundaries on the cavity length. Then the next 


assumed value of cavity length iss 


Sapper t Boer 
er ce 
Kit -— a 
The lower and upper boundaries on cavity length are 
determined by comparing (-%-), to (%) input. For (%) input 
greater than (ey the actual cavity length is greater than He, 


and so long as £; is greater than/f,, the lower boundary for the 


cavity length is then £; . For(%),,),; less than (%). 
£, vecomes the new upper boundary. Thus the boundaries on 
cavity length will converge to a solution. 

The following program listing uses this procedure. This 
program also uses the subroutines CPGEN, MATRIX and RMINV. 





96 
The iterative procedure for the super cavitating case 


was developed to overcome apparent nonconvergence of the method 
used in partial cavitation. However, the resulting nonconver- 
gence was based on theoretical cavity lengths only slightly 
greater than chord length (e.g << =~ 1.04). In this region 
of cavity length %/, increases very rapidly with cavity 
length. Convergence of any method for these cavity lengths 
will be slow. Therefore, the method used herein is probably 


an artificial result and not a required procedure, 
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Table C l 
Results of Cavity Length Computations 


_ — computed - No, of Geurst Percent 
cavity iterations ayy | difference 
length 


rm 
3 
3 
3 
6 
9 
6 
6 
hk 
8 
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